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ABSTRACT 


^  / 

■J 

This  paper,  which  is  the  first  in  a  two/part  study,  addresses  certain 
issues  concerning  the  small^train  theory  of  nonlinear  elasticity.  It  con¬ 
siders  isotropic  materials  which  possess  a  linear  response  in  shear  and  a 


nonlinear  response  in  dilatation,  and  (i)  establishes  an  explicit  necessary 
and  sufficient  condition  for  the  existence  of  piecewise  homogeneous  defor¬ 
mations,  (ii)  obtains  a  characterization  of  the  set  of  all  such  deforma¬ 
tions,  (iii)  derives  an  expression  for  the  ^^driving  traction^  on  a  surface 
of  discontinuity  in  the  strain,  and  finally  (iv)  discusses  the  notion  of  a 
kinetic  law.  While  the  analysis  is  carried  out  within  a  three-dimensional 


setting,  the  results  are  shown  to  have  a  particularly  simple  form  when 
expressed  in  terms  of  a  certain  constitutive  fxinction  .  In  Part  II  of 

this  study  we  examine  a  specific  boundary^value  prob^^OT/  ^  \ 
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1.  Introduction 

In  this  paper,  which  is  the  first  in  a  two-part  study,  we  show  that 
certain  features  of  the  finite  theory  of  elasticity  are  also  present  in  the 
small -strain  nonlinear  theory;  the  particular  class  of  constitutive  laws 
that  we  consider  here  is  one  that  has  been  used  to  model  the  mechanical 
response  of  ceramic  composites  undergoing  supercritical  phase  transforma¬ 
tions.  In  Part  II  we  will  examine  a  specific  boundary -value  problem. 

A  number  of  recent  studies  in  finite  deformation  elasticity  theory  have 
been  concerned  with  "nonelliptic  materials",  see  for  example  [1-10].  Such 
materials  are  capable  of  sustaining  deformations  whose  gradient  is  discon¬ 
tinuous  across  certain  surfaces  in  the  body;  this  leads  to  a  tremendous 
lack  of  uniqueness  of  solution  to  boundary-value  problems,  since  the  class 
of  functions  from  among  which  a  solution  is  sought  has  to  be  greatly 
enlarged  to  allow  for  such  deformations.  Horeover,  quasi- static  motions  of 
a  body  composed  of  such  a  nominally  elastic  material  can  involve  a  dissipa¬ 
tion  of  mechanical  energy  at  particles  located  on  a  moving  surface  of  dis¬ 
continuity,  [9] . 

Continuum  mechanical  treatments  of  stress -induced  phase  transforma¬ 
tions  in  solids  Involve  such  deformations,  e.g.  [6,7].  In  the  context  of 
phase  transformations,  a  surface  of  displacement  gradient  discontinuity 
corresponds  to  a  phase  boundary  separating  two  different  phases  of  the 
material,  and  the  aforementioned  non-uniqueness  might  be  thought  of  as 
arising  due  to  the  fact  that  the  classical  equations  of  the  continuum 
theory  do  not  account  for  the  kinetics  of  the  transformation. 
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In  the  present  study  we  examine  the  corresponding  issues  within  the 
infinitesimal  strain  theory  of  nonlinear  elasticity.  We  show  that  the 
aforementioned  phenomena  (of  discontlnuotis ,  dissipative,  non-unique  defor¬ 
mations)  persist  in  the  infinitesimal  strain  theory  too,  suggesting  that 
(in  some  sense)  it  is  the  constitutive  nonlinearity  rather  than  the  kine- 
matical  one  that  is  the  principal  source  of  these  features. 

In  this  study  we  restrict  attention  to  the  particular  class  of  consti¬ 
tutive  laws  that  were  proposed  by  Budiansky,  Hutchinson  and  Lambropou- 
los[14]  for  modeling  the  mechanical  response  of  certain  transforming  ceram¬ 
ics.  The  fracture  toughness  of  these  ceramic  composites  (which  contain 
second  phase  particles  that  undergo  a  phase  transformation)  was  known 
to  be  higher  than  that  of  the  brittle  ceramic  matrix  [11,12,13].  In  order 
to  model  this  phenomenon  at  the  continuum  level,  Budiansky,  Hutchinson  and 
Lambropoulos  [14]  derived  a  homogenized  constitutive  law  for  such  compos¬ 
ites  using  arguments  based  on  the  self-consistent  method.  They  argued  that 
since  the  transformation  leads  to  particles  twinned  into  layers  of  alter¬ 
nating  shear,  the  average  shear  associated  with  the  transformation,  from  a 
continuum  point  of  view,  is  essentially  zero.  Accordingly,  they  proposed 
(and  studied)  a  constitutive  law  with  a  linear  response  in  shear  and  a  tri- 
linear  response  in  dilatation;  see  also  Silling  [15] .  It  is  this  class  of 
materials  that  we  will  study  here  (modified  to  allow  the  dilatational 
response  to  be  arbitrary) . 

Chen  and  Reyes  Moral  [16]  have  experimentally  examined  the  relative 
importance  of  shear  and  dilatation  in  transforming  ceramics,  and  Lambrou- 
polos  [17]  has  proposed  a  more  general  constitutive  law  that  accounts  for 
both  of  these  effects.  Ve  do  not  consider  such  generalizations  here. 


In  this  paper,  we  first  recall  the  ellipticity  conditions  for  the 
three-dimensional  displacement  equations  of  equilibrium;  they  are  shown  to 
have  a  particularly  simple  interpretation  in  terms  of  the  stress  response 
function  of  the  material  in  uni-axial  deformation,  S(e).  Next,  in  Section 
3,  we  examine  conditions  under  which  a  three-dimensional  piecewise  homoge¬ 
neous  deformation  can  be  sustained  by  the  material,  and  derive  a  single 
necessary  and  sufficient  condition  for  the  existence  of  piecewise  homoge¬ 
neous  deformations.  This  condition  too  is  expressed  in  a  particularly 
simple  form  in  terms  of  Z(e):  in  addition  to  providing  infomation  on 
existence,  it  also  allows  us  to  characterize  the  set  of  all  possible  piece- 
wise  homogeneous  equilibrium  states. 

As  shown  by  Knowles [9],  when  the  theory  of  finite  elasticity  is  broad¬ 
ened  to  allow  for  equilibrium  fields  with  discontinuous  displacement  gra¬ 
dients  ,  the  usual  balance  between  the  rate  of  external  work  and  the  rate  of 
storage  of  elastic  energy  during  a  quasi-static  motion  no  longer  holds. 
Instead,  one  finds  that  mechanical  energy  may  be  dissipated  at  points  on 
the  surfaces  of  discontinuity.  This  in  tvim  permits  one  to  introduce  the 
notion  of  a  "driving  traction"  which  may  be  viewed  as  a  normal  traction 
that  the  body  applies  to  the  surface  of  discontinuity  at  each  of  its 
points.  In  Section  4  we  observe  that  a  dissipation  of  mechanical  energy 
can  also  occur  in  the  small -strain  theory  of  elasticity,  and  we  derive  an 
explicit  expression  for  the  driving  traction  in  the  case  of  the 
aforementioned  materials;  this  too  may  be  simply  expressed  in  terms  of 
2(e). 


The  stress  response  function  in  uni-axial  deformation  11(e)  plays  such' a 
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visible  role  in  all  of  these  results  because  (as  shown  in  Section  3)  the 
local  deformations  on  the  two  sides  of  a  surface  of  discontinuity  differ 
from  each  other  by  precisely  a  uni-axial  stretch  in  the  direction  normal  to 
that  surface. 

In  Section  5  we  briefly  discuss  the  need  for  additional  constitutive 
infomation  in  order  to  complete  the  theory.  As  discussed  there,  this 
might,  for  example,  take  the  form  of  a  "kinetic  law"  which  relates  the 
driving  traction  on  the  surface  of  discontinuity  to  its  velocity  of  propa¬ 
gation.  The  "flow  rule"  utilized  by  Budiansky  et  al[14]  is  equivalent  to  a 
particular  kinetic  law  as  will  be  discussed  more  fully  in  Part  II. 

The  results  in  this  paper  pertaining  to  the  existence  of  piecewise 
homogeneous  deforaations  (in  three -dimens ions)  have  a  similar  form  to 
analogous  results  for  isotropic.  Incompressible  elastic  materials 
undergoing  finite  plane  deformations,  [20].  Likewise,  the 
(three-dimensional)  driving  traction  formula  here  is  similar  to  the 
corresponding  formulae  for  finite  plane  and  anti -plane  deformations 
[20,21].  A  discussion  of  kinetic  relations  in  the  particular  setting  of 
the  one -dimensional  theory  of  bars  was  given  in  [22]. 

Finally,  in  Section  6  we  observe  that  the  driving  force  on  the  tip  of  a 
crack  is  generally  affected  by  the  presence  of  a  surface  of  strain  discon¬ 
tinuity,  even  if  the  crack  is  stationary;  see  [14].  A  relationship  between 
the  far  field  value  of  the  J- integral,  the  near- tip  value  of  J  and  the 
resultant  driving  force  on  the  surface  of  discontinuity  is  derived. 
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2.  Preliminaries. 

Consider  an  elastic  body  occupying  a  region  R  of  three-dimensional 
space.  Let  ^  be  the  position  vector  of  a  particle  in  R  and  let  u(x)  be  its 
displacement.  Suppose  that  there  Is  a  smooth  surface  S  which  lies  in  R, 
such  that  the  displacement  field  is  continuous  on  R  and  twice  continuously 
differentiable  on  R-S;  Vu  may  suffer  a  finite  jvunp  discontinuity  across  S. 
Let  H,  and  ^  denote  the  displacement  gradient  tensor,  the  infinitesimal 
strain  tensor  and  the  strain  deviator  respectively: 


HiJ 


'^i.j* 


l/2(  Ui  j  +  Uj^i). 

«ij  *  V3  ekk  ^ij> 


for  X  e  R-S. 


(2.1) 


Displacement  continuity  across  S  requires 


[[Uij]]ij  -  0  for  x,€  S  (2.2) 

for  all  vectors  ^  that  are  tangential  to  S  at  x;  [[•]]  indicates  the  jump 
across  the  surface  S.  Finally,  let  d(^  and  k(^  denote  the  respective 
strain  invariants  which  represent  the  dilatation  and  shear  at  a  particle 

A  -  tr  ] 

\  for  X  e  R-S.  (2.3) 

k  -  [2tr(e2)]V2,  J 


Next,  let  £(^  be  the  stress  tensor  field  on  R  and  suppose  that  £(x)  is 
continuously  differentiable  on  R-S;  ^  may  suffer  a  finite  jiuap  discontinu¬ 
ity  across  S.  Equilibrliim  in  the  absence  of  body  forces  requires 


‘"ij,j  -  0-  ‘^ij  -  ‘^ji 
[(ffij]]nj  - 


► 


(2.4) 


0 


for  €  R-S, 
for  jc  €  S, 


(2.5) 


where  n  is  a  iinlt  normal  vector  on  S.  A  surface  S  which  carries  jump  dis¬ 
continuities  in  and  ^  while  maintaining  displacement  and  traction 
continuity  is  called  an  eouilibrixim  shock  or  phase  boundary. 

Turning  to  the  constitutive  law  of  the  material  at  hand,  suppose  that 
it  is  homogeneous,  isotropic  and  hyperelastic.  The  elastic  potential  W 
then  depends  on  the  deformation  only  through  the  three  principal  invariants 
of  strain.  A  particular  case  of  special  interest  is  that  in  which  W 
depends  only  on  the  shear  and  dilatational  invariants  k  and  A: 

W(j^-W(k,A).  (2.6) 

The  stress-strain  relation  -  3V/d^  at  a  particle  x  €  R-S  then  specializes 
to 


-  (2A)aw/ak  +  (aw/aA  -  (2A/3k)aw/ak)  (2.7) 

If  the  material  is  such  that  the  mean  stress  (7£i/3  depends  on  the  deforma¬ 
tion  only  through  the  dilatation  one  can  show  using  (2.7)  that  it  is 

necessary  and  sufficient  that  (2.6)  have  the  separable  form  W(k,A)  -  f(k)  + 
g(A)  which  can  be  more  conveniently  written  as 

k  A 

W(k,A)  -  I  r(K)dK  +  I  a(?)d4  for  1^0,  -• «»<A<«;  r(0)-ff(0)-0.  (2.8) 

0  0 

(Alternatively,  one  can  show  that  the  components  of  deviatoric  stress 

depend  on  the  deformation  solely  through  the  components  of  deviatoric 
strain  if  and  only  if  W  has  the  form  (2.8).)  The  constitutive  functions 
r(k)  and  o(A)  may  be  readily  interpreted  as  follows:  in  a  simple  shear 


deformation  ui  •>  kxo,  uo  -  0,  U3  0,  the  shear  stress  component  012  i-s 
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found  from  (2.8),  (2.7)  to  be  a-^2  “  &  pure  dilatational  deforma¬ 
tion  Uj_  -  (A/3)x£,  one  finds  that  “  ff(A) .  Thus  the  function  r(k)  is 

the  shear _ stress  response  function  of  the  material  in  simple  shear,  while 

the  function  o(A)  is  the  mean  stress  response  function  of  the  material  in 
pure  dilatation. 

Finally,  we  further  specialize  (2.8)  to  the  case  in  which  the  shear 
stress  response  in  simple  shear  is  linear:  r(k)  -  nV..  This  is  motivated  by 
the  fact  that  such  constitutive  relations  appear  to  be  of  interest  in  the 
continuum  mechanical  modeling  of  certain  ceramic  composites  containing  par¬ 
ticles  which  undergo  stress  induced  phase  transformations,  (see  Budiansky 
et  al[14]).  Thus,  in  this  study  we  consider  materials  characterized  by  an 
elastic  potential 

A 

W(k,A)  -  (M/2)k2  +  I  for  k^O,  .•<:a«o,  (2.9) 

0 

where  n  (>0)  is  the  infinitesimal  shear  modulus  of  the  material.  The 
stress -strain  relation  (2.7)  now  specializes  to 

oij  -  2m  +  (a(A)  -  2MA/3)fiij.  (2.10) 

The  bulk  modulus  of  the  material  (2.9)  is 

B(A)  ■  o(A)/A  for  -«o<A<«o.  (2.11) 

It  is  useful  for  later  purposes  to  consider  the  response  of  this  body 
in  a  uni -axial  deformation  uj^  -  ex^,  U2  -  U3  -  0.  From  (2.10)  one  gets 
ail  “  where 


(2.12) 


The  displacement  equations  of  equilibrium  for  the  class  of  materials 
under  discussion  here  are,  by  (2.10),  (2.4),  (2.1),  (2.3) 

Cijkl(^^'ik,jl  -  0  for  R-S, 
where 

Oijkl^X/)  “  ^(^ik^lj  +  ^kj^li)  +  (^'(4)  -  2M/3)^lj^kl. 

The  system  of  partial  differential  equations  (2.13)  is 
(strongly)  elliptic  at  a  solution  u.  and  at  a  point  if 

Oijkl^t^j))  minjmijni  >  0  (2.15) 

for  all  unit  vectors  m  and  i^.  It  is  not  difficult  to  show  from  (2.14), 
(2.15)  and  /i  >  0  that  (strong)  ellipticity  prevails  if  and  only  if 

a'(A(j9)  >  -4/i/3  (2.16) 

where  A(^  -  (^(^(x)  is  the  dilatation  associated  with  the  given  deformation 
at  the  point  under  consideration.  Observe  from  (2.12)  that  this  elliptic* 
ity  condition  can  be  expressed  simply  in  terms  of  the  stress  response  func¬ 
tion  in  uni -axial  deformations  as 


(2.13) 


(2.14) 

said  to  be 


S' (A(x))  >  0.  (2.17) 

Thus  the  ellipticity  of  the  governing  equations  is  directly  related  to  the 
invertibility  of  the  stress  response  function  in  uni>axial  deformations. 
If  S  fails  to  be  monotonically  increasing  on  -«6  <  c  <  *,  ellipticity  will 
be  lost  at  some  deformation.  If  r'(c)  >  0  for  all  e,  we  say  that  the 
material  is  elliptic.  We  assume  throughout  that  S' (0)  >  0  so  that  ellip¬ 
ticity  prevails  at  the  undeformed  state;  since  S' (0)  -  >c+4/i/3  where  k  is 
the  infinitesimal  bulk  modulus,  this,  together  with  fi  >  0,  are  the  usual 
ellipticity  conditions  of  linear  elasticity. 


3.  Piecewise  homogeneous  displacement  fields. 

Not  all  homogenous,  isotropic  elastic  materials  characterized  by  the 
constitutive  relation  (2.10)  can  sustain  deformations  with  discontinuous 
strains.  In  this  section,  we  determine  a  simple  necessaxry  and  sufficient 
condition  on  the  material  which  determines  whether  or  not  it  can  sustain 
piecewise  homogeneous  deformations  of  this  type.  In  addition,  for  materi¬ 
als  that  can  sustain  such  deformations,  we  obtain  a  characterization  (in  a 
certain  sense)  of  the  entire  collection  of  possible  piecewise  homogeneous 
deformations . 

We  now  consider  the  special  case  in  which  R  coincides  with  all  of 

(x]^,X2,X3)- space,  S  is  a  plane  through  the  origin,  and  the  displacement 

gradient  is  constant  on  each  side  of  S.  Let  n  be  a  unit  vector  normal  to 

+  — 

the  plane  S,  and  let  R,  R  be  the  two  open  half- spaces  into  which  S  divides 
R  with  the  normal  pointing  into  R.  The  field  equations  (2.13)  will  then 
be  trivially  satisfied  in  R  -S,  and  all  that  remains  to  be  fulfilled  are 


the  Jiimp  conditions  (2.2)  and  (2.5). 


Consider  the  piecewise  homogeneous  displacement  field 


u  -  •< 


+  + 

H  X  for  x€R  , 


(3.1) 


I  «  S  for  j^eR  . 

+ 

where  the  displacement  gradient  tensors  ^  and  ^  are  constant  and  distinct: 

(3.2) 


+ 

H  >«  H. 


+  —  — 

Define  c ,  e ,  A  and  A  by 


«ij  -  1/2  (H^j  +  Hji) 

+  + 

^  “  «kk* 


‘Ij  -  V2  (Hij  +  Hji) 

^  ”  «kk- 


(3.3) 

(3.4) 


The  displacement  field  (3.1)  will  be  continuous  across  S  If  and  only  If 

+ 

ij  -  fj  for  all  unit  vectors  ^normal  to  n,  (3.5) 

while  by  (2.5),  (2.10)  the  tractions  will  be  continuous  across  S  If  and 
only  if 

+  +  +  -  -  - 
2m  ‘ij^j  +  (<r(A)  -  2MA/3)n£  -  2m  ‘Ij^j  +  (a(A)  -  2MA/3)ni.  (3.6) 

+ 

Given  a  tensor  ,  the  shock  problem  consists  of  finding  a  tensor  JH  and  a 
unit  vector  n  such  that  (3.5)  and  (3.6)  (with  (3. 3), (3. 4))  hold. 

Ue  first  establish  a  necessary  condition  which  must  hold  If  the  shock 

+ 

problem  Is  to  have  a  solution.  To  this  end,  suppose  that  given  H,  there  Is 


a  tensor  H  and  a  imit  vector  ii  such  that  (3. 3) -(3. 6)  hold.  It  can  be 
readily  shown  that  (3.5)  holds  if  and  only  If  there  exists  a  vector  &  such 
that 

+ 

^ij  "  ^ij  ^i’^j  •  (3.7) 

Thus  (3.3),  (3.4),  (3.7)  yield 
+ 

fij  “  +  l/2(ainj  +  ajnj^) ,  (3.8) 

+ 

A  -  A  +  aif^i-  (3-9) 

Turning  next  to  the  requirement  (3.6)  and  multiplying  it  by  the  components 
of  any  unit  vector  normal  to  n  gives 

+ 

eij^inj  -  ‘ij-^i^j*  (3.10) 

which  in  view  of  (3.8)  simplifies  to 

a^ii  -  0.  (3.11) 

Since  this  must  hold  for  all  unit  vectors  in  the  plane  S,  it  follows 
that  a  is  parallel  to  xi: 

a  -  on.  (3.12) 

By  (3.12),  (3.9), 

+ 

a  -  A  -  A.  (3.13) 

Moreover  (3.8)  can  be  written,  in  view  of  (3.12),  as 
+ 

«ij  “  «ij  +  “"i"j  • 


(3.14) 


Finally,  multiply  the  traction  continuity  condition  (3.6)  by  n^  and  use 
(3.14),  (3.13)  to  obtain 

+  +  -  - 

o(A)  +  4/iA/3  -  a(A)  +  4/iA/3 ,  (3.15) 

which,  in  terms  of  the  uni-axial  deformation  response  function  Z(A) ,  reads 
+ 

2(A)  -  2(A) .  (3.16) 

Next  we  will  show  that  if  the  (necessary)  condition  (3.16)  holds,  then 

this  in  fact  guarantees  the  existence  of  a  solution  to  the  shock  problem. 

+  + 

In  order  to  show  this,  suppose  Chat  ii  is  a  given  tensor.  Define  A  by 

•  + 

(3.4)]^,  (3.3)^.  If  there  exists  a  number  A  (f*  A)  such  that  (3.16)  holds, 

then  (for  each  arbitrary  unit  vector  n)  we  can  define  a  by  (3.13),  a, by 

-  + 
(3.12)  and  by  (3.7).  It  may  be  readily  verified  that  these  tensors  H  and 

H  automatically  satisfy  the  requirements  (3.5),  (3.6)  of  displacement  and 

traction  continuity  across  Che  plane  with  unit  normal  n.  Thus  we  have  the 

following  result: 

Proposition:  Given  a  tensor  E,  there  exists  an  associated  piecewise 
homogeneous  equilibrium  shock  if  and  only  if  there  is  a 
number  A  (p<  ^  ■  ^kk^  such  that  (3.16)  holds. 

When  the  constitutive  law  is  such  that  the  stress  response  function 
2(e)  is  monotonically  increasing  (in  %diich  case  the  material  is  elliptic) 
we  see  from  the  preceding  proposition  that  the  material  cannot  sustain  a 
piecewise  homogeneous  deformation.  On  the  other  hand,  if  the  material  is 
such  that  2' (A)  0  on  some  interval,  then  since  2' (0)  >0,  it  follows  chat 

piecewise  homogeneous  deformations  will  exist  for  suitably  chosen  values  of 

+ 

H. 


In  adiitlon  to  providing  InformaClon  concerning  the  existence  of  a 
piecewise  homogeneous  deformation  associated  with  the  given  displacement 
gradient  H,  the  preceding  result  also  peralts  us  to  characterize  the  set  of 
all  such  deformations  which  can  be  associated  with  that  H:  it  states  that 
for  every  number  A  for  which  (3.16)  holds,  and  for  all  choices  of  the  unit 
normal  vector  n,  one  can  construct  an  acceptable  H.  Let  E  denote  the 
following  set  in  the  (A,  A) -plane: 

+  -  +  -  +  • 

S  -  {  (A.  A)  I  S(A)  -  S(A),  A  A  ).  (3.17) 


According  to  the  preceding  proposition,  given  a  displacement  gradient  ten- 

sor  H,  the  associated  shock  problem  has  a  solution  if  and  only  if  there  is 

-  +  -  +  - 
a  number  A  such  that  (A,  A)  €  E  where  A  -  moreover,  all  tensors  li 

+  - 
that  can  be  connected  to  H  by  a  shock  are  generated  by  all  numbers  A  for 

•t 


which  (A,  A)  e  S.  The  set  S  characterizes  the  collection  of  all  possible 
shocks.  A  sketch  of  the  curve  S  in  the  (A,  A) -plane,  corresponding  to  a 
particular  class  of  materials^ will  be  given  in  Section  5. 


Finally,  we  note  that  according  to  (3.7),  (3.12),  (3.13)  the  displace- 

+■  - 

ment  gradient  tensors  and  H  are  related  by 
•  +  .  ^ 

"  Hjj  +  (A  -A)  n^nj  ;  (3.18) 

this  implies  that  the  deformation  on  R  is  equivalent  to  the  deformation  on 
+ 

R  together  with  a  uni-axial  stretch  in  the  direction  normal  to  the  shock 
surface.  This  is  presumably  the  reason  ^y  the  stress  response  function  in 
uni-axial  deformation  S(e)  plays  such  a  central  role  in  the  preceding  (and 
subsequent)  results. 


4.  Driving  Traction 

We  now  consider  a  quiasl- static  motion  of  the  body  and  let  u(*,t),  tQ  < 
t  <  t]^,  be  a  one-parameter  family  of  solutions  of  the  displacement  equa¬ 
tions  of  equilibrium  (2.13)  of  the  type  described  in  Section  2.  Let  S^-  C  R 
be  the  family  of  shocks  associated  with  this  motion,  and  assume  that  the 
particle  velocity  v(x,t)  -  au(x,t)/at  exists  and  is  continuous  in  (x, t)  for 
x^  e  R-S^.  tg  <  t  <  tj^,  and  that  v  is  piecewise  continuous  on  R  x  [tg.tj^]. 

Let  d(t)  denote  the  difference  between  the  rate  of  external  work  (on 
any  fixed  regular  region  II  c  R)  and  the  rate  at  which  elastic  energy  is 
being  stored  (in  H) : 


<l(t) 


J  ‘^Ij^'j^i 

an 


dA 


d_ 

dt 


I  »(£)  dV, 

n 


to  <  t  <  ti; 


(4.1) 


d(t)  is  the  rate  of  dissipation  of  mechanical  energy  in  the  region  H.  By 
adapting  to  the  present  small -strain  theory  the  analysis  given 
by  Knowles  [9] ,  one  can  show  that  d(t)  may  be  written  as 

d(t)  -  f  n«V  dA,  (4.2) 

J 

s^nii 


where  f(x, t)  is  defined  by 


f  -  •  n,  for  x  €  Sf  tQ  S  t  ti,  (4.3) 


P(x, t)  is  the  energy-momentum  tensor 


Pij  “  *  ‘^kj  R'St.  to  <  t  <  t]i  ,  (4.4) 

and  V(x, t)  is  the  velocity  of  a  point  on  the  moving  surface  S^.  If  the 
motion  happens  to  be  smooth,  jP(*,t)  will  be  continuous  across  and  so 
(4.2),  (4.3)  gives  d(t)  -  0  for  tg  <  t  <  tj^.  .  In  general  however  the 
dissipation  rate  d(t)  f*  0  whenever  n  intersects  S^.. 


If  we  postulate  that  at  each  instant,  the  rate  of  storage  of  energy  in 
n  cannot  exceed  the  rate  of  external  work  on  n,  then  we  must  require  the 
dissipation  rate  d(t)  to  be  non-negative  for  all  sub-regions  n  and  all 
instants  t.  Thus,  by  (4.2), 

f  Vjj  >  0  for  xeS^.,  (4.6) 
where  is  the  normal  velocity  of  a  point  on  the  surface  : 


~  Z’  Z  to<t<ti.  (4.7) 

Alternatively,  the  dissipation  Ineoualitv  (4.6)  can  be  shown  to  be  a 
consequence  of  the  second  law  of  thennod3mamlcs  under  Isothermal  condi¬ 
tions;  see  [9].  In  general,  given  an  equilibrivm  state,  the  inequality 
(4.6)  restricts  the  direction  in  which  the  surface  may  move  in  a  quasi- 
static  motion  commencing  from  this  state. 

A  particularly  simple  expression  for  the  driving  traction  can  be 
derived  in  the  case  of  materials  characterized  by  the  special  elastic 
potential  (2.9).  First,  from  (3.7),  (3.12)  one  has 
+ 

[[Hki]]  -  -  aninij.  (4.8) 

Next,  in  view  of  (2.10),  (2.12),  (4.8)  and  the  continuity  of  traction, 

+  +  +  + 
n  “  -  o(  2M«ijn£nj  -  IftA  +  S(A)  ).  (4.9) 

However  from  (2.3),  (2.1)  and  (3.14)  one  obtains 
+  +  +  - 

[[k2]]  -  -4a«ijn£nj  -  2o2  +  2o(A  +  A)/3,  (4.10) 

+ 

which  can  be  used  to  eliminate  the  term  e^jn^nj  from  (4.9)  to  give 

+  +  -  -  +  + 

[[‘^kjHki«inj]l  -  -  k2)/2  +  2m(a2  -  a2)/3  -  02(4).  (4.11) 

Finally,  since  (2.9)  and  (2.12)  provide 

A 

W(<^  -  /ik2/2  -  2mV3  +  I  S(?)  d^, 

0 


(4.12) 


L8-- 


I 

equations  (4.3),  (4.4),  (4.11),  (4.12)  and  (3.13)  yield  Che  desired  expres¬ 
sion 


+ 

4 

f  +  +  . 

f  -  I  S(4)  dA  -  S(A)  (A  -  A)  for  JceS^, 

A 


(4.13) 


for  Che  driving  traction.  In  the  finite  theory,  formulae  of  this  general 
form  have  been  derived  In  the  special  case  of  "normal  shocks"  In  plane  and 
anti-plane  finite  deformations  of  Isotropic,  Incompressible  elastic  solids, 
[19,20]. 

It  Is  useful  to  write  (4.13)  as 
+  - 

f  -  F(A,  A)  for  t(ySt:Sti  ,  (4.14) 

where  F  Is  the  fvinctlon  defined  on  the  set  S  by 

+ 

A 

F(A,  A)  -  I  Z(A)  dA  -  2(A) (A  -  A)  for  (A,  A)6  2.  (4.15) 

A 

By  (4.14),  the  driving  traction  f  at  a  point  on  the  phase  boundary 

depends  only  on  the  local  dilatations  A,  A  on  the  two  sides  of  S^;  f  does 

-t*  - 

not  depend  on  Che  amounts  of  shear  k,  k,  nor  on  the  orientation  of  S^. 
Moreover,  In  view  of  (4.15)  and  (3.16),  the  value  of  f  may  be  Interpreted 
geometrically  as  the  difference  between  the  area  xmder  the  uni -axial 
deformation  stress-strain  curve  between  A  and  and  the  area  of  the 

rectangle  on  the  same  base  with  height  2(X) . 


5.  Kinetic  Law.  An  example. 

In  Part  II,  we  will  present  an  example  which  shows  that  boundary -value 
problems  formulated  in  the  conventional  manner,  for  materials  characterized 
by  (2.9),  may  suffer  from  a  tremendous  lack  of  uniqueness.  This  is  known 
to  be  the  case  in  the  finite  theory  as  well  (e.g.  [18,19]).  This  non-uni¬ 
queness  suggests  that  the  theory,  as  formulated,  is  deficient,  and  that  it 
ought  to  be  supplemented  with  additional  constitutive  information.  One v way 
in  which  to  implement  this  is  to  postulate  a  constitutive  relation,  or 
"kinetic  law",  which  applies  to  particles  on  S^,  and  relates  the  driving 
traction  f  to  the  normal  velocity  of  propagation  Vn  of  the  phase  boundary. 

In  order  to  formalize  this,  let  f){  and  fg^  be  the  supremum  and  infemum 
of  the  function  F(A,  A)  on  the  set  E.  Then,  one  might  suppose  that  there 
is  a  constitutive  function  V(*)  defined  on  [fg,,  f^j]  such  that 

Vn  -  V(f)  on  St  .  tQ  S  t  S  ti.  (5.1) 

In  order  to  conform  to  the  dissipativlty  inequality  (4.7),  V  must  be  such 
that 

f  V(f)  >0  for  f  €  [fg,,  f^].  (5.2) 

The  form  (5.1)  is,  of  course,  merely  an  example  of  a  class  of  kinetic  laws 

that  might  be  imposed;  it  could  be  generalized  to  include  dependence  on 

other  local  variables  as  well  so  that,  for  example,  a  kinetic  law  might 
+  - 

read  -  V(A,  A),  where  the  constitutive  function  V  is  defined  on  S. 

In  order  to  illustrate  this  (and  some  of  the  preceding  results)  we  now 
choose  the  dllatational  stress  response  function  in  the  constitutive  law 
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(2.10)  to  be  as  follows: 


^(A) 


for  0  <  A  <  A^, 

•  ^A  +  aj(A-Aj4)/(Am-Aj{)  for  A^  <  A  <  Ag,, 

.  ^A  +  oj  for  A  >  Ag,; 


(5.3) 


Aq,  Ajf  and  a<j>  are  material  constants  such  that 


/3>0,  AgP»AM>0, 

/8dg,  +  a-j;  >  0, 

<  ~ot- 


(5.4) 


The  second  condition  In  (5.4)  Implies  that  ff(Agj)  >  0,  while  (5.4)3  ensures 
that  the  system  of  equations  (2.13)  Is  non-elllptlc  when  A^  <  A(x)  <  Ajj. 
In  this  example  we  will  confine  attention  to  the  range  A  >  0  and 
consequently  we  have  left  o(A)  iindeflned  for  negative  values  of  Its 
argviment.  The  specific  constitutive  law  (5.3),  (5.4)  Is  the  one  considered 
by  Budlansky,  Hutchinson  and  Lambropoulos  [14]  in  the  case  of  super¬ 
critical  transformations.  As  shown  In  Figure  1,  as  the  dilatation 

Increases,  the  mean  stress  first  rises  linearly  to  a  SMxlmum  value  - 

^A){,  It  then  declines  linearly  to  the  val\ie  9g|  -  ^Ag  -f  oj,  and  finally 
rises  again  with  the  Initial  slope 


"The  response  function  of  this  material  In  uni -axial  deformations  Is 
given  by  (5.3),  (2.12)  as 


S(0 


ae  for  0  ^  c  ^  A)(, 

-  ac  +  9j(c-AK)/(Ag|-^)  for  A){  £  c  ^  Ag, 

06  +  Oj  for  6  't  A|g, 


(5.5) 


where  we  have  set 


a  -  fi  +4m/3. 


(5.6) 


Finally,  we  Introduce  the  following  additional  notation  which  pertains 
to  certain  special  points  on  the  stress-strain  curve  shown  in  Figure  1: 


^1  "  ^ 
^ol  “ 


-  An  -  £7^/0, 

^o3  “ 


.1 


(5.7) 


Note  that  the  straight  lines  which  join  (^M3 ^  • 
(Aoi,a(Aoi))  to  (Ao3,a(Ao3)) ,  and  (Amj^,o(Anji) )  to  (A^.a^),  each  have  the 
same  slope  -4/i/3;  see  (3.15)  for  the  significance  of  this.  Moreover,  AqJ^ 
and  Aq3  are  seen  to  obey  the  conditions 


where 


S(Aoi)  -  Z(Ao3)  -  (i:„+%)/2, 

“  2(Am)  ,  ZIn  •  S(An)  . 


1 


(5.8) 


The  set  E  for  this  material,  which  characterizes  the  complete  set  of 
possible  shocks  in  the  (A,  A) -plane,  may  be  readily  found  from  (3.17), 
(5.5).  It  consists  of  the  points  fiQ  the  polygon  ABCDEFA  shown  in  Figure  2, 
except  for  the  vertices  A  and  D  which  lie  on  the  line  A  -  A. 

We  turn  next  to  equation  (4.15)  which  defines  the  driving  traction 
function  F  on  this  set  S.  Explicit  formulae  for  F  may  be  readily  derived 
from  (4.15),  (5.5).  For  example,  when  (A, A)  €  EF  ,  one  finds 

F(A,  A)  -  (-ot/«){  -  (Zb+%)/2).  (5.9) 

We  do  not  display  the  remaining  formulae  here.  It  Is  particularly  useful 
to  know  the  sign  of  the  driving  traction,  since  then  the  direction  of  pro- 
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pagatlon  of  the  phase  boxindary  is  known  through  the  dissipatlvity 
inequality  (4.6).  The  sign  of  F  may  be  read  off  from  (5.9)  (and  the 
analogous  formulae  appropriate  to  the  other  points  on  S) ;  one  finds  that 

>  0  for  (A.  A)  e  (AB]+[BP)+(DE]+[EQ), 

F(A,  A)  -  >  <0  for  (t,  A)  €  (PC]+[CD)+(QF]+[FA) ,  (5.10) 

-  0  for  (A,  A)  -  P  or  Q. 

(The  s}niibol  (AB]  in  (5.10)  denotes  the  set  of  all  points  on  the  line  AB 
excluding  the  end  point  A  but  including  the  point  B.)  The  points  P  and  Q 
which  are  associated  with  zero  driving  traction  are  sometimes  refered  to  as 
"Maxwell  states".  They  are  given  by 

+■  - 

(A.  A)  -  (Ao3.  Aoi)  and  (Aoi.  ^o3'>  •  (5.11) 

where  A^]^  and  A03  were  defined  in  (5.7);  see  also  (5.8).  Also,  one  finds 

that  the  driving  traction  achieves  its  largest  value  f^j  at  B  (and  also  at 

E)  and  its  smallest  value  f,|  at  C  (and  also  at  F) .  These  values  are 

fM  -  -  Za)/2a  (>0) ,  (5.12) 

f m  -  ax(%  -  2:„)/2a  «0) .  (5.13) 

Finally,  Figure  3  shows  an  exasiple  of  a  kinetic  function  V  that  might 
be  used  in  the  kinetic  law  (5.1).  It  is  consistent  with  the  admissibility 
requirement  (5.2).  In  the  example  which  will  be  discussed  in  Part  II,  we 

will  see  how  in  a  specific  problem,  the  kinetic  relation,  together  with  an 

initiation  criterion,  can  be  \ised  to  resolve  the  non -uniqueness  referred  to 
earlier.  In  that  example  we  will  find  that  the  kinetic  relation  of  Figure 
3  generally  leads  to  rate -dependent  "viscoplasticity- like"  response.  Two 
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special  cases  which  lead  to  reversible,  dissipation- free  response  and  to 
rate - independent  plasticity- like  respone  will  also  be  discussed  there. 


6.  Concluding  remark:  driving  force  on  a  crack- tip 

In  this  section  we  briefly  comment  on  the  driving  force  on  a  crack- tip 
when  the  crack  is  contained  in  a  body  composed  of  the  material  (2.10).  For 
simplicity,  suppose  that  the  body  is  a  slab  containing  a  traction-free 
through-crack  (Figure  4)  and  that  the  loading  is  such  that  the  deformation 
is  planar.  Suppose  further  that  the  body  is  composed  of  the  material 
(2.10)  with  the  constitutive  function  S(<)  defined  by  (2.12)  being  non¬ 
monotone.  By  the  analysis  in  Section  3,  deformations  of  this  body  can 
Involve  shocks.  Suppose  for  definiteness  that  there  is  a  single  (cylindri¬ 
cal)  phase  boundary  S  as  sho«m  in  Figure  4;  C  Is  the  curve  along  which  S 
Intersects  the  (X]^,X2) -plane.  The  deformation  is  smooth  at  all  points  of 
the  body  inside  C  (excluding  points  on  the  crack  itself)  as  well  as  at  all 
points  outside  C. 


Let  Fq  and  be  two  closed  curves  as  shown  In  Figure  4  with  Fg  being 
entirely  within  C  and  F*  entirely  outside.  The  values  of  the  J- Integral 
associated  with  these  two  curves  are  respectively, 


'^tip 


To 


ds- 


(6.1) 


where  ds  denotes  arc  length,  jn  Is  the  unit  outward  normal  vector  on  the 
appropriate  curve,  and  are  the  components  of  the  energy-momentum 
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tensor : 


W  ^0)8  -  <^7^  ^70* 


(6.2) 


(Greek  subscripts  take  the  values  1  and  2  only.) 


The  J- integral  is  path- independent  nrovi 


C;  this,  together  with  the  traction- free 


nature  of  the  crack  surface  yields  the  alternate  expression 


-  f 


Pl^  ds, 


(6.3) 


where  are  the  limiting  values  of  as  a  point  on  C  is  approached  from 
within.  Similarly 


- 


■t* 

Pl^  ds, 


(6.4) 


Combining  (6.3)  and  (6.4)  gives 


dtip  “ 


(6.5) 


Next,  in  view  of  (6.2),  (2.1),  displacement  continuity  (2.2),  and 

traction  continuity  (2.5),  one  sees  that 


[(Po^]]_n;8  <*o  “  0  o"  C, 


(6.6) 


where  j^is  a  unit  tangent  vector  on  C.  Thus  the  vector  [[P,,^]]  n^  is  normal 


--25-- 

to  the  curve  C,  and  by  (4.3), 

-  f  on  C;  (6.7) 

f  is  the  driving  traction  on  the  shock.  Finally,  combining  (6.5)  with 
(6.7)  provides  the  desired  expression 

“^tip  ”  '^00  f  f  n]^  ds.  (6.8) 

C 

Equation  (6.8)  states  that  the  driving  force  Jtip  on  the  crack- tip 
equals  the  difference  between  J*  (the  "applied  value  of  J")  and  the 
resultant  driving  force  on  the  shock.  Thus  in  general,  Jtip  J®-  (This 
was  also  noted  by  Silllng[15] . )  In  certain  exceptional  cases,  for  example 
if  the  deformation  is  such  that  f  -  constant  on  C,  the  integral  in  (6.8) 
will  vanish  and  then  Jtip  “  J«-  The  value  of  the  shock  driving  traction  f 
depends  on  (and  is  determined  by)  the  particular  kinetic  relation  governing 
the  evolution  of  the  shock.  If  the  resultant  driving  force  on  the  shock  is 
in  the  positive  Xi-direction  then,  by  (6.8),  Jpip  <  J®- 
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Figure  2.  The  set  E  characterizing  all  possible  shock  states. 


Figure  4.  Geometry  of  cracked  slab  with  phase  boundary  C  and 
Integration  paths  r^,  . 
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